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ABSTRACT 

Using the number and sizes of observed gravitational lenses, I derive upper limits on the dark matter 
content of elliptical galaxies. On average, dark matter can account for no more than 33% of the total mass 
within one effective radius {Re) of elliptical galaxies, or 40% of the mass within 2Re (95% confidence 
upper limits). I show that galaxies built from Cold Dark Matter (CDM) mass distributions are too 
concentrated to comfortably satisfy these limits; a high-density {Hm — 1) CDM cosmology is ruled out 
at better than 95% confidence, while a low-density, flat cosmology is only marginally consistent with the 
lens data. Thus, lensing adds to the evidence from spiral galaxy dynamics that CDM mass distributions 
are too concentrated on kiloparsec scales to agree with real galaxies, and extends the argument to 
elliptical galaxies. 

Lensing also provides a unique probe of the very inner regions of galaxies, because images are predicted 
to form near the centers of lens galaxies but are not observed. The lack of central images in deep maps of 
radio lenses places strong lower limits on the central densities of galaxies. The central densities of CDM 
galaxies are too low on '--^lO parsec scales. Supermassive black holes can help suppress central images, 
but they must lie well off the observed black hole-bulge mass correlation in order to satisfy current limits 
on central images. Self-interacting dark matter, or any other modification to regular cold dark matter, 
must simultaneously reduce the densities on kiloparsec scales and increase the densities on parsec scales 
in order to satisfy the unique constraints from lensing. 



Oh; 
O 



X 



1. INTRODUCTION 

Cuspy mass distributions are a robust prediction of the 
popular Cold Dark Matter (CDM) paradigm. Numerical 
simulations of coUisionless cold dark matter predict den- 
sity profiles with p oc r~" and a ~ 1.0-1.5 at small radii 
(e.g., Navarro, Frenk & White 1996, 1997; Moore et al. 
1998, 1999; Jing & Suto 2000; Klypin et al. 2000), and 
this prediction does not depend on particular cosmogo- 
nies or initial conditions (Huss, Jain & Steinmetz 1999a, 
1999b) or on the specific form of the dark matter power 
spectrum (Eke, Navarro & Steinmetz 2000). Adding dis- 
sipative baryons makes mass distributions even more con- 
centrated (e.g., Blumenthal et al. 1986; Dubinski 1994). It 
is important to compare the predicted mass distributions 
with real galaxies to test the CDM paradigm as the ex- 
planation for the formation and growth of structure in the 
universe. 

The dynamics of spiral galaxies have provided the most 
extensive observational tests of CDM. Salucci (2001; also 
see Salucci & Burkert 2000, and references therein) sug- 
gests that normal spiral galaxies have dark matter ha- 
los with large constant-density cores as opposed to cusps. 
Debattista & Sellwood (1998) and Weiner, Sellwood & 
Williams (2001) argue that fast-rotating bars require 
dark matter densities lower than predicted by CDM. The 
slowly-rising rotation curves of dwarf galaxies and low 
surface brightness galaxies also seem to imply constant- 
density cores (e.g., Flores & Primack 1994; Moore 1994; 
McGaugh & de Blok 1998; Blais-Ouellette, Amram & 
Carignan 2001; de Blok et al. 2001). However, several au- 
thors have argued that the Hi rotation curves may have ap- 
peared artificially shallow due to beam smearing, and that 



better data are actually consistent with cuspy CDM mass 
distributions (van den Bosch et al. 2000; van den Bosch 
& Swaters 2000; Swaters, Madore & TrewheUa 2000). In 
rebuttal, McGaugh & de Blok (1998) and de Blok et al. 
(2001) claim that beam smearing does not affect their con- 
clusions; and Moore (2001) argues that only two of the 19 
galaxies analyzed by van den Bosch & Swaters (2000) are 
actually consistent with CDM. To summarize, many argue 
that CDM mass distributions are too concentrated to agree 
with the dynamics of spiral galaxies, but the conclusions 
are still subject to some vigorous debate. 

Other tests of CDM mass distributions should not be af- 
fected by beam smearing. Navarro & Steinmetz (2000) and 
Eke et al. (2000) consider the global dynamical properties 
of spiral galaxies in terms of the Tully-Fisher (TF) relation 
between luminosity and circular velocity. They find that 
CDM can reproduce the slope and scatter of the TF rela- 
tion, but has some trouble with the zero point. In a high- 
density (Om = 1) cosmology, CDM model galaxies are too 
concentrated to agree with the TF zero point, while in a 
low-density cosmology the models are marginally consis- 
tent with the data. Taking an entirely different approach, 
Rix et al. (1997) study the line-of-sight velocity profile 
of the elliptical galaxy NGC 2434 using detailed dynami- 
cal models. They find that the galaxy is consistent with 
CDM models, but the strength of the conclusion is limited 
by systematic uncertainties such as the orbital anisotropy. 
These tests do not indicate a fundamental problem with 
CDM mass distributions, but they do not strongly favor 
the CDM models, either. 

The dynamical tests are fundamentally limited by the 
need to interpret data from luminosity distributions before 
drawing conclusions about mass distributions. Given the 
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importance of the CDM paradigm in modern cosmology, 
it is desirable to develop additional tests that are inde- 
pendent of, and hopefully less ambiguous than, the dy- 
namical tests. One excellent possibility is gravitational 
lensing, because it offers a direct probe of mass distribu- 
tions. Individual gravitational lenses robustly determine 
the masses of individual galaxies, and the statistical prop- 
erties of the lens sample constrain properties of the galaxy 
population (e.g., Maoz & Rix 1993; Kochanek 1993, 1995, 
1996; Cohn et al. 2001). Individual lenses and lens statis- 
tics both imply that elliptical galaxies (which dominate 
the mass-selected sample of lens galaxies; e.g., Kochanek 
et al. 2000) have approximately isothermal mass distri- 
butions out to several kiloparsecs, in agreement with the 
evidence from dynamics (e.g., Rix et al. 1997) and X-ray 
elliptical galaxies (e.g., Fabbiano 1989). 

It is not clear whether the mass distributions implied 
by lensing are consistent with the predictions of CDM. 
Most lensing studies, even those that consider a wide range 
of density profiles (e.g., Kochanek 1995; Barkana 1998; 
Chae, Khersonsky & Turnshek 1998; Cohn et al. 2001; 
Muhoz, Kochanek & Keeton 2001), consider only single- 
component mass models. However, lens galaxies are likely 
to have at least two components (a stellar galaxy and 
a dark matter halo), which may contribute comparable 
amounts of mass (e.g., Rix et al. 1997). Both compo- 
nents are necessary to explain the distribution of lensed 
image separations and the fact that galaxies arc much bet- 
ter lenses than more massive groups of galaxies (Keeton 
1998; Porciani & Madau 2000; Kochanek & White 2001). 
Only by allowing two components can we use lensing to di- 
rectly test whether real galaxies are consistent with CDM 
mass distributions. 

The goal of this paper is to use two-component 
star-|-halo models for lensing to test the CDM paradigm. 
The focus is on elliptical galaxies because they dominate 
lens statistics. The outline of the paper is as follows. Sec- 
tion 2 defines the models, and Section 3 reviews the lensing 
calculations. In Section 4, the observed number and sizes 
of lenses are used to evaluate the global properties of the 
models. In Section 5, lensing is used to examine the very 
inner regions of galaxies. Finally, Section 6 offers a discus- 
sion and conclusions. 

2. STAR-I-HALO MODELS 

This section defines star-l-halo models for elliptical 
galaxies in the context of the CDM paradigm. Section 2.1 
discusses models for the stellar and dark matter compo- 
nents, and Section 2.2 gives normalizations for the mod- 
els. Only spherical models are considered, because they 
are sufficient for calculations of the number and sizes of 
lenses. Dcpartiircs from spherical symmetry mainly af- 
fect the relative numbers of 2-image and 4-image lenses 
(see Keeton, Kochanek & Seljak 1997; Rusin & Tegmark 
2001). 

2.1. Model components 

A simple model for the stellar components of elliptical 
galaxies is the Hernquist (1990) model, which has a density 
profile 



This profile is described by a scale radius r,, but the pro- 
jected profiles of elliptical galaxies are usually described 
by an effective (or half- mass) radius i?e; they are related 
by Vg = 0.551i?e- The total mass of a Hernquist model is 
M ~ The projected surface mass density, in units 

of the critical density for lensing, is 

m _ nm _ ^ (2 + x^).F(a;)-3 
R/vg, Kg = PsTs/^cri and the function T{x) 



where x 
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The critical density is Ecr = {c^Ds)/{^'KGDiDia) where 
Di and Dg are angular diameter distances to the lens and 
source, respectively, and Dis is the angular diameter dis- 
tance from the lens to the source (e.g., Schneider, Ehlers & 
Falco 1992). The gravitational deflection for a Hernquist 
model is (see eq. 14 below) 



bR{R) = 2Ksr, 



x[l-T{x) 
a;2 - 1 



(4) 



Navarro, Prenk & White (1996, 1997, hereafter NFW) 

have argued that dark matter halos foTind in cosmologi- 
cal A^-body simulation of coUisionless dark matter have a 
"universal" density profile of the form 

Ps 



p{r) 



(5) 



(r/r,)(l + r/r,)2 ' 

where Vs is a scale radius and ps is a characteristic density. 

It is convenient to replace the scale radius with a "concen- 
tration" parameter C = r2oo/rs, where r2oo is the radius 
within which the mean density of the halo is 200 times 
the critical density of the universe, which is often taken to 
mark the boundary of a relaxed halo (e.g., Crone, Evrard 
& Richstone 1994; Cole & Lacey 1996; Navarro et al. 1996, 
1997). The characteristic density is then 

_ 200 

Ps-^ Pcrit(z) in(l + C)- C/{1 + C) ' ^ ' 

where Pcrit{z) is the critical density of the universe at the 
redshift of the halo. The lensing properties of an NFW 
model are given by Bartelmann (1996). 

More recently, Moore et al. (1998, 1999; also see Jing & 
Suto 2000; Klypin et al. 2000) have argued that the cen- 
tral regions of simulated halos are steeper than the NFW 
profile. They advocate a density of the form 

Ps 



p{r) 



(7) 



(r/r«)i-5[l + (r/r,)i-5] ' 
For Moore halos, I define a concentration parameter C = 
'"20o/?'(-2) in terms of the radius ?'(-2) at which the loga- 
rithmic slope of the density is —2. This definition is equiv- 
alent to the definition of the concentration for NFW halos, 
and Keeton & Madau (2001; also see Wyithe, Turner & 
Spcrgcl 2000) argue that it is the best generalization of 



the concentration. The radius 



r(-2) 



is related to the scale 



radius Vg in eq. (7) by t"(_2) = 0.630r.,. The characteristic 
density of a Moore halo is 



Ps = 25pcritiz) 



ln[l-FC3/2/2] 



(8) 



CDM AND STRONG LENSING 



o 



> 

CD 
O 



2 
1.5 
1 

0.5 



2 




T — I — I — I — I — I — I — I — I — I — I — I — I — r 



C = 5. fooo.=0-15 



^C=10, U,=0.15 - 



total 



halo 

Initial NFW 



I I I I I I I I I I I I 

-C = 5, U,=0.10 



I I I I I I I I I I I 




-C=15, feooi=0.15 - 



I I I I I I I I I I I I I I 

-C=15. fcooi=0.10 - 



radius r/R. 



Fic;. 1. Rotation curves for sample star+halo models with NFW halos. Each panel has the specified values of the concentration C of the 
initial halo and the cooled mass fraction fcool\ models have -Re/''200 = 0.03. The solid curves show the total rotation curves, while the 
dotted and dashed curves show the contributions from the galaxy and halo, respectively. For comparison, the long dashed curves show the 
rotation curves of the initial NFW halos before adiabatic contraction. The velocities are scaled by the peaJs velocity of the galaxy component. 



The dark matter halo models were derived from studies 
of coUisionless dark matter. They do not hold in the pres- 
ence of dissipative baryons, because as the baryons cool 
and condense into a galaxy they modify the gravitational 
potential and thus the dark matter distribution (e.g., Blu- 
menthal et al. 1986; Dubinski 1994). Fortunately, there is 
a simple analytic prescription called adiabatic contraction 
for computing the changes to the dark matter distribu- 
tion; it seems to agree well with gasdynamical simulations 
(e.g., Blumenthal et al. 1986; Flores et al. 1993), even 
for merger scenarios thought to produce elliptical galaxies 
(Gottbrath 2000). Appendix A gives an analytic solution 
for adiabatic contraction of an arbitrary dark matter halo 
by a Hernquist galaxy. Adiabatic contraction depends on 
the mass ratio of the cooled galaxy component to the to- 
tal virial mass, fcooi = Mgai/Mtot, which is presumably 
no larger than the global baryon fraction of the system, 
fbar = Mbar/Mtot = ^b/^M- (There may be baryons 
that remain hot and distributed throughout the halo, so 
fcool < fbar-) The virial mass and radius of the system 
factor out to provide overall scalings, so the solution also 



depends on the fraction fcool of the concentration C of the 
initial halo, and the effective radius of the galaxy (specifi- 
cally Re/r2oo)- 

To illustrate the star-l-halo models, Figure 1 shows ro- 
tation curves for various values of the parameters. For a 
fixed stellar component, decreasing fcool increases the to- 
tal mass of the halo {Mtot = fcooi ^ ^gai), which raises 
the rotation curve. Increasing the concentration of the 
initial halo packs more of the dark matter into the inner 
regions of the system, which also raises the inner rotation 
curve. In other words, changing either parameter affects 
the amount of mass contained within a few effective radii 
of the galaxy. Lensing can distinguish between the two 
parameters only if it is sensitive to the detailed shape of 
the galaxy mass profile inside a few Rg. 

Figure 1 offers two important qualitative results. First, 
the galaxy and halo components can easily combine to pro- 
duce a rotation curve that is relatively flat from ~ 0.5Re 
to several Re- In other words, star+halo models can nat- 
urally produce net mass distributions that are fairly close 
to p a throughout much of the galaxy. Second, com- 
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paring the rotation curves of the halo before and after 
adiabatic contraction illustrates that the modification by 
the baryons can significantly increase the halo mass within 
a few Re, especially for less concentrated halos. Figure 2 
also shows that adiabatic contraction affects NFW profiles 
more dramatically than Moore profiles, especially for large 
fcooU which tends to reduce the differences between NFW 
and Moore model galaxies. 
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radius that has somewhat larger scatter but is easier to 
use. Schade, Barricntos & Lopez-Cruz (1997) find that at 

= the M{B)-logRe relation is 

MAB{B)-5logh = -3.331og(i?eA"^kpc)-18.15±0.06, 

(9) 

and the relation evolves with redshift in a way that is con- 
sistent with the fading of stellar populations due to passive 
evolution. Early-type dwarf galaxies, on the other hand, 
appear to form a population that is disjoint from giant 
galaxies. Binggeli & Cameron (1991, 1993) demonstrate 
this effect in the Virgo cluster, and a fit to their data yields 

Mab{B) = -13.1 log(i?e/kpc) - 13.7 , (10) 
although with significant scatter. The break between giant 
and dwarf galaxies occurs somewhere around an absolute 
magnitude of —16 or —18, but it is not sharp. The ex- 
act location of the break has little effect on lens statistics 
because these low-mass galaxies contribute little lensing 
optical depth. 

Simulated dark matter halos do not all have the same 
profile. Halos of a given mass have a range of concentra- 
tions; and cluster-mass halos are systematically less con- 
centrated than galaxy- mass halos (e.g., Jing & Suto 2000; 
Bullock et al. 2001). Bullock et al. (2001) characterize the 
scatter by fitting NFW profiles -'^ to simulated halos and ob- 
taining a set of concentration parameters consistent with 
the log-normal distribution 

1 r {logic /Cmed{M,z)]f 



p{logC\M,z) 



radius r/R 



exp< - 



2a2 
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Fig. 2. — A comparison of rotation curves for NFW and Moore 
models. The four panels show the various components of the rota- 
tion curve. The solid curves indicate NFW models and the dot- 
ted curves show Moore models. Results are shown for C = 5, 
fcool = 0.15, and Re/r200 = 0.03. 

2.2. Normalizations 

The CN0C2 field galaxy redshift survey (Lin et al. 1999, 
2001) gives the luminosity function of early-type galaxies 
at redshifts 0.12 < z < 0.55. The luminosity function 
is parametrized as an evolving Schechter (1976) function, 
and Lin et al. (2001) give parameter values for two cos- 
mologies: a high-density, flat universe with matter density 
Qm = 1; and a low-density, flat universe with matter den- 
sity r^M = 0.2 and cosmological constant ^\ = 0.8. I use 
a Hubble constant Hq = 50 km s~^ Mpc~^ for the ^Im = 1 
cosmology to mimic the Standard Cold Dark Matter cos- 
mology, and Ho = 65 km s~^ Mpc~^ for the = 0.2 flat k{R) 
cosmology. I convert the luminosity of the stellar com- 
ponent into a mass using population synthesis models by 

Bruzual & Chariot (1993), modehng early-type galaxies ~ 2i? 

with an old coeval stellar population. 

I place a dark matter halo around each galaxy and use 
empirical correlations to normalize the galaxies and ha- 
los. Bright early-type galaxies are observed to populate a <j)ii{R) 
"fundamental plane" in the space of surface brightness, ef- 
fective radius, and velocity dispersion, with very little scat- 
ter away from this plane (e.g., Djorgovski & Davis 1987; ~ ^ 
Dressier et al. 1987). Projecting out the velocity disper- P 
sion yields a correlation between luminosity and effective fn 

^ Bullock et al. (2001) remark that it would be possible to fit other profiles to halos but argue that the eq. (11) captures the full range of halo 
properties seen in their simulations. 
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where gq = 0.18, and the median concentration varies 
systematically with mass and redshift as Cmed{M,z) oc 
M~^/^(l -|- z)~^. The scatter in halo properties is impor- 
tant for lensing because more concentrated halos arc much 
better lenses. To include this effect, I use halos drawn ran- 
domly from eq. (11). I normalize the distribution in terms 
of the parameter C defined to be the median concentration 
of 10^2 /j-i Mq halos at redshift z = 0. The value of C is 
predicted by simulations (see §4.2), but I take it to be a 
free model parameter. 

3. LENSING METHODS 

The adiabatic contraction solution gives the mass profile 
M(r) of the final system. The system's projected surface 

density k{R) = Y.{R)/Ticr and lensing deflection (^ij can 
then be written as (see Kceton 2001) 

m'(r) 



K200 



^200 



dr 



dy 



1 



1 + ?r 



X 



K200 ''200 R 
1^200 T 200 



JR 



to' [r^/i 

m{r) 



r^^Jr^ — R? 



,2^3/2 



(12) 
(13) 

(14) 
(15) 
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where the radii arc written in units of r2oo, the mass 
m(r) = M(r)/Af2oo is written in units of the total mass 
inside r2oo, and m'(r) = dm{r)/dr. Eqs. (13) and (15) 
represent variable transformations that give the integrals 
a finite range, which is useful for numerical integration. 
The strength of the system as a gravitational lens is mea- 
sured by the dimensionless parameter 

K200 = 9 ^ (16) 



7rr; 



200 • 



0.00467 



M. 



200 



-,1/3 


\H{z)] 




[ Ho \ 



-,4/3 



IQio/i-iMeJ [ Ho J rnD, ' 

where rn = c/Hq is the Hubble distance. This parame- 
ter is the mean projected surface density of the system in 
units of the critical density for lensing. In general, K200 
is considerably less than unity because most halos can act 
as strong gravitational lenses only in a high-density region 
near the core, not all the way out to the virial radius. 

The images corresponding to a given source are found 
by solving the lens equation, 

u = R- MR) , (17) 
whore u is the angular position of the source relative to 
the lens (see Schneider et al. 1992 for a full discussion). 
The magnification of an image at position R is 

t,{R) = {l-MR)-\l-<l>RR)-K (18) 

Here (j)RB = d{(j)R)/dR, which can be computed efficiently 
using the identity <pR + 4>RR = 2k. In general, a spher- 
ical lens has two radii at which the magnification is infi- 
nite. These radii correspond to "critical curves" in the 
image plane, which map to "caustics" in the source plane. 
The outer or tangential critical curve lies at the Einstein 
ring radius of the lens; a source directly behind the lens 
produces a ring image with radius Ve- The inner or radial 
critical curve lies at a small radius Vcr- The source posi- 
tion corresponding to an image at Vcr, which I label Ugut, 
marks the boundary of the region where lensing yields mul- 
tiple images. (The equations for rcr and Uout are given in 
Appendix B.) A source with u < Uout has three images, 
one outside te, one between ve and rcr, and one inside 
rcr', the innermost image is usually demagnified and unde- 
tected (see §5). A source with u > Uout has a single image, 
which is outside r^. 

Computing the statistics of gravitational lenses requires 
summing over populations of lenses and sources, and ac- 
counting for "magnification bias," or the fact that a flux- 
limited survey may include lenses where the source is in- 
trinsically fainter than the flux limit, but lensing magnifi- 
cation brings the object into the sample (e.g.. Turner 1980; 
Turner, Ostriker & Gott 1984). The number of lenses with 
a total flux greater than S expected to be found in a survey 
with lensing selection functions described by T is 



where Zs is the source redshift, dV is the comoving vol- 
ume element (see, e.g., Carroll, Press & Turner 1992), and 
dn/dM is the mass function of halos that can serve as 
lenses. The integral over C incorporates the scatter in 
halo properties defined in eq. (11). The factor J^{u) indi- 
cates whether a lens associated with a source at u would 
be detected given the selection functions. The volume, 
mass, and concentration integrals sum over the popula- 
tion of possible lens galaxies. The Zs integral allows for a 
distribution of source redshifts, where [dNsrc{> S)/dzs] dzg 
is the number of sources brighter than flux S that lie in the 
redshift range Zg to Zg + dzg- Finally, the distribution of 
image separations is found by computing dN\ens/dO, and 
mean quantities are found by averaging over the predicted 
lens population. 



4. THE NUMBER AND SIZES OF LENSES 

In this section, the global properties of the star-|-halo 
models are evaluated using two quantities from lens statis- 
tics: the number of lenses, or more specifically the fraction 
of sources that are multiply imaged; and the distribution of 
lensed image separations.^ Section 4.1 reviews the data. 
Section 4.2 presents results for a fiducial set of models, 
while Section 4.3 considers systematic effects including the 
source redshift distribution, the galaxy formation redshift, 
and the density profile. 

4.1. Data 

More than 50 galaxy-mass lenses are known, and their 
properties have been compiled by the Cf A/ Arizona Space 
Telescope Lens Survey (CASTLES; see Kochanek et al. 
2001). This sample includes lenses from a variety of sur- 
veys as well as serendipitous discoveries; thus, the par- 
ent (or source) population is unknown, and the CASTLES 
sample cannot be used to test the number of lenses. By 
contrast, the distribution of image separations in the sam- 
ple probably can be used, because it is insensitive to the 
size of the source population. 

The largest homogeneous statistical survey for lenses 
is the Cosmic Lens All-Sky Survey (CLASS; Helbig 
2000; Browne 2001). The sample comprises 10,499 flat- 
spectrum radio sources with flux S' > 30 mJy at 5 GHz, 
and the flux distribution can be described as a power law 
dNsrc/dS oc S'^ with ly w —2.1 (see Rusin & Tegmark 
2001). The survey includes 18 lenses, all of which have 
image separations 6 < 3", and the survey is believed to 
be complete at image separations 0.3" < 6 < 15" (Helbig 
2000; Phillips et al. 2001). Because this paper focuses on 
lensing by elliptical galaxies, I omit two CLASS lenses that 
are known to be produced by spiral galaxies (B 0218-1-357 
and B 1600+434). For the CLASS lenses where the lens 
galaxy type is not known, I assume an elliptical galaxy 
because most lens galaxies are ellipticals (e.g., Kochanek 
et al. 2000) and because this is the conservative approach 
(as shown below). In the CLASS survey both the source 
and lens populations are known, so the sample can be used 
to test both the number of lenses and the distribution of 
image separations. 

^ A third interesting quantity is the ratio of four-image lenses to two-image lenses, which can be used to constrain the angular shape of lensing 
mass distributions (e.g., Kochanek 1996; Rusin & Tegmark 2001). This test requires non-spherical lens models, and it is not very sensitive to 
the mass profile of lensing halos. 
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Fig. 3. — Image separation histograms for the CLASS data (solid lines) and for sample models (dotted lines). Model results are shown for 
the fiducial models in an Q,m = 0.2 flat cosmology. The model parameters are indicated in each panel. 



The number of lenses can be tested (the "N test" ) by us- 
ing Poisson statistics to compare the CLASS sample with 
predictions from the models. The image separations can 
be tested (the test") by using a Kolmogorov-Smirnov, 
or K-S, test (e.g., Press et al. 1992) to compare the ob- 
served and predicted distributions of image separations. 
The CLASS and CASTLES samples can both be used for 
the 9 test, although the tests are not independent because 
the CASTLES sample contains the entire CLASS sample. 
The two samples are consistent in the sense that a K-S 
test does not reveal a significant difference between their 
separation distributions. The CASTLES sample is larger 
and thus less sensitive to statistical peculiarities (such as 
the lack of CLASS lenses with 9 > 3", perhaps). The 
CLASS sample, however, has better information about 
source fluxes and redshifts. In an attempt to compromise 
between the two samples, I perform the 9 test with both 
samples and conservatively adopt the weaker of the two 
results. 

The largest uncertainty in the models arises from the 
source redshift distribution, which is not known for the 
full CLASS sample. Marlow et al. (2000) report redshifts 
for a small subsample of 27 sources from the CLASS sam- 
ple. They find a mean redshift of (z) = 1.27, which is 
comparable to that found in other radio surveys at compa- 
rable fluxes (Drinkwater et al. 1997; Henstock et al. 1997; 
Falco, Kochanek & Mufioz 1998). They also find evidence 
for a difi^erence between the galaxy and quasar popula- 
tions in the sample, with (zgai) = 0.18 for 8 galaxies and 
(zqso) = 1-72 for 19 quasars. It is not clear at this point 
whether the subsample fairly represents the full sample. 

^ Recall that the calculation explicitly includes scatter in the halo 
concentration. 



To examine possible systematic effects, I consider a set 
of models with all sources placed at the mean redshift of 
the subsample, and an alternate set of models with source 
redshifts distributed according to the subsample. 

4.2. Basic results 

Consider a fiducial set of models in which the halos be- 
fore adiabatic contraction are modeled with NEW profiles, 
the galaxies have old stellar populations that formed at 
redshift Zf = 5, and all the sources are placed at the mean 
redshift of the CLASS spectroscopic subsample, Zg = 1.25. 
Figure 3 compares model predictions with the data from 
the CLASS sample. As the median-^ concentration C in- 
creases or the cooled mass fraction fcooi decreases, the 
number of lenses increases and the distribution of image 
separations shifts to higher values. Physically, increasing 
C or decreasing fcooi raises the amount of dark matter in 
the inner parts of halos (see §2.1), leading directly to more 
and larger lenses. 

Comparing the models to the data using the A'^ and 9 
statistical tests yields confidence limits on the model pa- 
rameters, as shown in Figure 4. There is a band in the 
upper left of the {C, fcooi) plane where the models are 
consistent with both the number of lenses and the distri- 
bution of image separations. Moving to larger C or smaller 
fcooi increases the number and sizes of predicted lenses. In 
the hatched region, the predicted lenses are generally too 
big, and the models are ruled out (at 95% confidence) by 
the test. In the cross-hatched region, the models are fur- 
ther excluded because they predict too many lenses (the 

properties (see eq. 11), so the models are charax;terized by the median 
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Flc. 4. Confidence regions in the (C,fcool) plane for the fiducial models. The shaded regions below the diagonal curves are excluded 
at 95% confidence by the lens data; the lower and upper curves correspond to the A'^ and 9 tests, respectively. The shaded regions above 
the horizontal lines are excluded by measurements of the cosmic baryon density H;,. The lower curve corresponds to = 0.019 ± 0.0024 

from measurements of deuterium (Tytler et al. 2000), and the upper curve corresponds to ^li,h? < 0.037 (95% confidence) from the cosmic 
microwave background (Tegmark et al. 2001). The arrows on the x-ax.es indicate concentrations predicted by CDM simulations (see text). 
Results are shown for two cosmologies. 



A'' test).^ The contours from the N test are based on the 
assumption that 16 CLASS lenses are produced by clHpti- 
cal galaxies. If the number of CLASS lenses with elliptical 
galaxies turns out to be smaller than 16, the N contours 
would move further up and to the left, strengthening the 
constraints from lensing. 

There is little difference between the lensing constraints 
in the two cosmologies shown in Figure 4, which seems sur- 
prising because it is traditionally argued that lens statis- 
tics are quite sensitive to a cosmological constant A (e.g., 
Turner 1990; Kochanek 1996). The traditional argument 
is based on models where the lens galaxy population is 
obtained by taking the local comoving number density of 
galaxies and assuming that it holds out to redshift z 1 
in all cosmologies; in this case, the number of lenses is 
very sensitive to the volume of the universe to z ~ 1, and 
hence to A. By contrast, my models arc based on counts 
of galaxies at z ~ 0.5 from the CN0C2 field galaxy red- 
shift survey (Lin et al. 1999, 2001). In these models, the 
volume factor required to convert from number counts to 
number density (or luminosity function) essentially can- 
cels the volume factor that appears in the lensing analysis; 
the number of lenses is roughly proportional to the num- 
ber counts of galaxies, and is not very sensitive to A. In 
other words, using models normalized by number counts 
of galaxies at z ~ 0.5 makes lens statistics only weakly 
sensitive to cosmology. 

CDM simulations make specific predictions about the 



concentration: in a cosmology with Qm = 1, (? — 11 

for standard CDM, or (7 ~ 7 for tilted CDM where the 
power spectrum has shape parameter F = 0.2; and in an 
0.M = 0.2 flat cosmology, (7^8 (Navarro et al. 1997). 
These values are indicated by arrows on the x-axes in Fig- 
ure 4. Cosmic baryon censuses give limits on fcooi- Be- 
cause Jcooi gives the fraction of a system's mass that has 
cooled into the baryonic galaxy, it is a lower limit on the 
baryonic content of the system and hence should not ex- 
ceed the cosmic baryon fraction, Vth/Q.M (e.g.. White et 
al. 1993). The upper limits on fcooi derived from mea- 
surements of using the cosmic microwave background 
(e.g., Tegmark, Zaldarriaga & Hamilton 2001) and the 
deuterium/hydrogen ratio and big bang nucleosynthesis 
(e.g., Tytler et al. 2000) are indicated by horizontal lines 
in Figure 4. Some elliptical galaxies contain hot. X-ray 
emitting gas that is probably primordial gas that never 
cooled; the cool stellar component may contain as little 
as half of the baryons (e.g., Brighenti & Mathews 1998). 
The presence of hot gas would reduce the upper limit on 
fcooi to something below fib/f^M; but because the actual 
amount of gas and its presence across the galaxy popula- 
tion are not well understood, I focus on the conservative 
upper limit from fJb. 

The lens data reject models where concentrated, mas- 
sive dark matter halos make elliptical galaxies overly eSi- 
cient lenses — a large portion of the {C, fcooi) plane. At 



At small C and large fcooh the models predict too few lenses that are too small compajred with the data. These constraints apply beyond 
the upper left corner of the (C, fcooi) plane in Figure 4. 
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the concentrations found in CDM simulations, Icnsing re- 
quires baryon fractions that are incompatible with a high- 
density universe. The lens data are formally compatible 
with a low-density universe, but only in a narrow corner 
of parameter space where galaxy-mass halos must be very 
efficient at cooling their baryons. The general conclusion, 
then, is that galaxies constructed from CDM mass distri- 
butions are too concentrated to agree with lens statistics, 
especially in a high-density CDM cosmology. 

There is clearly a degeneracy between C and fcooi in 
Figure 4, which is not surprising because both parame- 
ters affect the central mass that determines the lensing 
properties. It is therefore interesting to define an integral 
quantity, 

which is the ratio of halo mass to galaxy mass inside some 
radius r, where the average is over the lens population. 
Note that M is defined using the mass in spheres. The 
mass ratio allows model-independent statements about the 
dark matter distribution in early- type galaxies, which are 
summarized in Table 1. In the fiducial models, dark mat- 
ter can account for up to 29-33% of the mass inside Re 
(95% confidence upper limit), and up to 35 40% of the 
total mass inside 2Re. In other words, dark matter can 
contribute a moderate fraction of the mass in the inner re- 
gions of elliptical galaxies, but it is not the dominant mass 
component at small radii. These interesting upper limits 
result from the distribution of image separations; observed 
lenses are too small to be consistent with larger dark mat- 
ter contributions. The lensing limits are similar to but 
stronger than those derived from a dynamical analysis of 
the nearby eUiptical galaxy NGC 2434 (Rix et al. 1997). 
They arc consistent with the lower limits on dark matter 
in ellipticals derived from the relationship between X-ray 
temperature and stellar velocity dispersion (Loewenstein 
& White 1999). 







Table 1 






IIalo/Galaxv Mass 


R ATK) 




Radius 


n&i = 1 


Om = 0.2 flat 


Case 1 


Re 


M < 0.50 


M < 0.41 




2Re 


M < 0.66 


M < 0.55 


Case 2 


Re 


M < 0.43 


M < 0.27 




2Re 


M < 0.57 


M < 0.35 


Case 3 


Re 


M < 0.50 


M < 0.44 




2Re 


M < 0.68 


M < 0.60 


Case 4 


Re 


M < 0.50 


M < 0.41 




2Re 


M < 0.65 


M < 0.52 



Note. — 95% confidence upper limits on the halo/galaxy 
mass ratio, defined in eq. (20), computed at two radii for 
two cosmologies. The four different cases are defined in 
the text. 

4.3. Systematic effects 

There are three systematic effects that may be impor- 
tant for the models. Figure 5 shows how the results from 
the fiducial models (case 1, Figure 5a) are changed by each 
effect, and Table 1 gives the updated constraints on the 
halo/galaxy mass ratio A4. In case 2 (Figure 5b), the 



fixed source redshift is replaced by a redshift distribution 
to match the CLASS spectroscopic subsample (Marlow et 
al. 2000). The new models predict larger image separa- 
tions, so the region excluded by the 6 test stretches up 
and to the left. While the entire (C, fcooi) plane is now 
formally excluded either by lensing or by fib, it is not clear 
how strongly to interpret this result, because the CLASS 
spectroscopic subsample may not fairly represent the full 
CLASS sample. The strength of the conclusions will ul- 
timately be limited by the extent to which the redshift 
distribution of the full CLASS sample can be determined. 
Nevertheless, it is important to discover that the redshift 
distribution may actually worsen the discrepancy between 
the data and the models. 

In case 3 (Figure 5c), the redshift at which the stel- 
lar populations formed is reduced from Zf = f> to Zf = 3. 
The yomigcr stellar populations have smaller mass-to-light 
ratios, so the galaxies are poorer lenses (because the lu- 
minosity function is held fixed). Thus, the models pre- 
dict fewer and smaller lenses, and the regions excluded by 
lensing move down and to the right in the ((7, fcooi) plane. 
Nevertheless, the changes in the lensing constraints are 
small; lens statistics are not very sensitive to the galaxy 
formation redshift, provided that the stellar populations 
of elliptical galaxies are old. The Fundamental Plane of 
elliptical galaxies in rich clusters out to z = 0.83 (e.g., van 
Dokkum et al. 1998) and of elliptical lens galaxies in low- 
density environments out to z ~ 1 (Kochanek et al. 2000) 
indeed implies old stellar populations, Zf > 2. 

Finally, in case 4 (Figure 5d), the initial NFW profiles 
are replaced by steeper Moore profiles. Moore halos have 
more mass in the central regions than NFW halos, even 
for a fixed concentration parameter, and thus yield bet- 
ter lenses. Hence, the models predict more and larger 
lenses, and the excluded regions move up and to the left 
in the {C, fcooi) plane. The change is not very dramatic, 
however, because of the effects of adiabatic contraction. 
Moore halos, which are denser than NFW halos to be- 
gin with, experience a smaller density enhancement under 
adiabatic contraction (see Figure 2). In other words, adi- 
abatic contraction tends to erase some of the differences 
between NFW and Moore models. 

These results suggest that systematic effects do not 
weaken the discrepancy between models and data, and 
may even strengthen it. CDM star-|-halo models are at 
best marginally consistent with the statistics of strong 
lenses, and may be quite inconsistent depending on the 
distribution of source redshifts in the full CLASS sample. 
As a relatively model-independent conclusion, the lensed 
image separations imply that dark matter can contribute 
no more than about 33% of the total mass inside Re, or 
about 40% of the mass inside 2Re (95% confidence; see 
Table 1). CDM halos appear to be too concentrated to 
agree comfortably with this constraint. 

5. ODD IMAGES 

In this section, the very inner regions of star-|-halo mod- 
els are evaluated with lensing. The fact that most lenses 
do not show the expected central or "odd" images places 
strong lower limits on the central densities of galaxies. Sec- 
tion 5.1 reviews the data. Section 5.2 presents results, and 
Section 5.3 offers a discussion. 
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Fig. 5. — Confidence regions in the (Cj/coof) plane, shown for the four cases defined in the text, in an Qm = 0.2 fiat cosmology. Panel (a) 
is the same as Figure 4b but is included here for completeness. 



5.1. Data 

It can be proved mathematically that a single thin lens 
with a smooth (i.e., non-singular) projected mass density 
and a finite mass always producics an odd number of im- 
ages (Burke 1981; also see Schneider et al. 1992). In other 
words, lenses are generally expected to have three or five 
images, but are usually observed to have two or four. The 
apparent paradox is resolved by noting that for lenses with 
high central densities, one of the images is close to the cen- 
ter of the lens and demagnified.^ If the central density is 



high enough, the central image may be highly dcmagni- 
fied and therefore very difficult to detect. For example, 
Appendix B shows that for an isothermal sphere with a 
small core radius, or for a power law density p oc r~" with 
a « 2, the mean magnification of central or "odd" images 
can be quite small. Every two- or four-image lens may 
therefore be a three- or five-image system where one of 
the images remains undetected. 

Odd images are expected to be rare in optical observa- 
tions, because they would be swamped by light from the 
lens galaxies. The only lens with an odd number of optical 



^ Alternatively, if the projected mass density is singular the odd image theorem formally breaks down. Odd images still appear, though, 
provided the central density cusp is shallower than p oc (see Appendix B). 
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Fig. 6. Contours of the fraction of (two-image) lenses where the odd image is brighter than 1% of the brightest image (fodd ^ 0.01). 
Results are shown for an Qm = 0.2 flat cosmology. The left panel shows models with NFW halos (case 1), and the right panel show models 
with Moore halos (case 4). 



images is APM 08279+5255 (Ibata et al. 1999): the third 
image is either a standard odd image, in which case it re- 
quires a shallow density cusp a < 0.4 for p oc r~" (Mufioz 
et al. 2001), or else it represents a special image configu- 
ration produced by an edge-on disk (Keeton & Kochanek 
1998). Radio observations should be much more sensitive 
to odd images, because the lens galaxies (as opposed to 
the sources) are rarely radio loud. The only candidate odd 
image detected in the radio is in MG 1131+0456 (Chen & 
Hewitt 1993), although the possibility that the lens galaxy 
is radio loud cannot be ruled out in this case. 

The CLASS lens sample ofi'ers high-resolution and high- 
dynamic range (noise level ^50 /iJy beam^^) radio maps 
of lenses with compact radio sources, and thus should be 
quite sensitive to odd images. Consequently, the fact that 
no CLASS lens shows an odd image leads to strong up- 
per limits on how bright the odd images can be. Rusin 
& Ma (2001) tabulate the upper limits on six two-image 
CLASS lenses. To factor out the unknown brightness of 
the source, they quote upper limits in terms of the flux ra- 
tio fodd defined to be the fiux of the odd image relative to 
the fiux of the brightest image. The 5cr upper limits range 
from fodd < 0.0083 for B 0739+366 (Marlow et al. 2001) to 
fodd < 0.00049 for B 0218+357 (Biggs et al. 1999). Nor- 
bury et al. (2001) give hmits on odd images for CLASS 
lenses with more than two images. When quantifying odd 
images relative to other lensed images, four-image lenses 
are much more sensitive to asymmetry in the lens galaxy 
than two-image lenses. Hence, I restrict attention to the 
two-image CLASS lenses where spherical models are suffi- 
cient for interpreting odd images. 



5.2. Results 

For power law models with p oc r^^ ° the mean magnifi- 
cation of odd images is unity, so odd images are not highly 
dcmagnified (sec eq. B7 in Appendix B). This simple pre- 
diction does not strictly apply to star+halo models, be- 
cause even in the initial halos the broken power laws aflfect 
lensing via projection, and adiabatic contraction increases 
the central density. However, it does suggest that the odd 
images in star+halo models arc worth investigating. The 
best way to draw conclusions from observational limits on 
odd images is to use models of individual systems that 
take into account not only the detection limits but also 
constraints on the global lens model from the observed 
images (e.g., Cohn et al. 2001; Munoz et al. 2001; Rusin 
& Ma 2001; Norbury et al. 2001). Instead, in this sta- 
tistical analysis I examine the distribution of odd images 
predicted by star+halo models to understand the general 
trends (c.f., Wallington & Narayan 1993). 

The upper limits on odd images in the CLASS dou- 
bles range from fodd < 0.0083 down to fodd < 0.0005; to 
be very conservative, we can simply say that no CLASS 
double has an odd image brighter than fodd = 0.01. In 
contrast. Figure 6 shows that star+halo models predict 
that more than 20% of lenses should have odd images with 
fodd > 0.01; and for the range of parameters allowed by 
the number and sizes of lenses, the fraction is more like 
30%. In other words, star+halo models predict that de- 
tectable odd images should be quite common, in confiict 
with observations. This result is not terribly sensitive to 
the model parameters. It is surprisingly similar for NFW 
and Moore models, despite the differences in the dark mat- 
ter cusps. The explanation is again adiabatic contraction: 
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Fig. 7. — The cumulative fraction of (two-image) lenses where the ratio of the odd image to the brightest image is greater than /odd- 
Results are shown for models with C = 7.7 and fcool = 0.19, in an Hm = 0.2 flat cosmology. The galaxy components are modeled as 
p oc r~"{rs + r)™"*, and each curve shows results for a particular value of a; the fiducial Hernquist model corresponds to a = 1.0. The 
initial mass distribution is modeled with an NFW (left panel) or Moore (right panel) profile. The heavy dashed curves show the upper limits 
derived from six two-image CLASS lenses (Rusin & Ma 2001). 



most of the mass in the cores of the galaxies was puhed in 
by adiabatic contraction, which affects NFW halos more 
strongly than Moore halos and thus tends to reduce the 
differences between the two models. The fraction of lenses 
with detectable odd images is a strong function of the de- 
tection threshold (see Figures 7 and 8), so the discrepancy 
between data and models grows as the upper limit on fodd 
is decreased. 

Odd images arc sensitive to the density profile at small 
radii, so the assumption of a Hernquist model galaxy with 
a p oc cusp should be examined. Faber et al. (1997), 
Ravindranath et al. (2001), and Rest et al. (2001) find 
that the surface brightness distributions I{R) of early-type 
galaxies have a range of cusps; luminous early- type galax- 
ies have cores or shallow cusps (/ oc with 7 < 0.3, 
corresponding to p oc r~" with a < 1.3), while fainter 
galaxies have steeper power-law cusps. Hence, the Hern- 
quist model seems reasonable for the massive galaxies that 
dominate lensing. Still, for completeness 1 consider models 
where the Hernquist galaxy is generalized to an arbitrary 
cusp using the density profile p oc r~"{rs+r)°'~'^ ■ Figure 7 
shows the results for models with C = 7.7 and fcooi = 0.19 
(a point close to the 9 boundary in Figure 4). Steep cusps 
suppress odd images, but only if they are considerably 
steeper than the luminosity cusps in luminous early-type 
galaxies. As a corollary, steep cusps also make galaxies 
more efficient lenses and thus aggravate the discrepancy 
between the models and the observed number and sizes 



of lenses. In other words, cusps do not provide a very 
attractive resolution to the odd image problem. These 
conclusions apply to both NFW and Moore dark matter 
models. 

The simple star-|-halo models may not be sufficient for 
this analysis, because many galaxies are observed to con- 
tain central supermassive black holes (e.g., Magorrian et 
al. 1998; Gebhardt et al. 2000a, 2000b; Ferrarese & Merritt 
2000; Merritt & Ferrarese 2001). While black holes have 
little effect on the number and sizes of lenses (they barely 
affect the potential on kiloparsec scales), they may sup- 
press or even eliminate odd images (Mao, Witt & Koop- 
mans 2001). To understand their effects on lens statis- 
tics, I add black holes to the star+halo models, where 
the mass of the black hole is determined from the veloc- 
ity dispersion of the galaxy using the empirical correlation 
Mbh = (1.3O±O.36)xlO^M0 (CT/200km s-i)4-72±o.36 (pgj._ 
rarese & Merritt 2000; Merritt & Ferrarese 2001). Fig- 
ure 8 shows that black holes normalized by this relation 
have little effect on odd images down to fodd 0.001. 
Making the black holes systematically more massive in- 
creases the suppression, but only for the faintest odd im- 
ages. Black holes must lie off the Ferrarese & Merritt re- 
lation by at least a factor of 10 in mass before they begin 
to affect odd images at the fodd = 0.01 level. These con- 
clusions again apply with almost equal strength to NFW 
and Moore models. 



^ The correlation was derived for nearby galaxies, but for simplicity I use it at all redshifts. This approach is conservative if black holes grow 
no faster than their surrounding galaxies, as in the models by Haehnelt & Kauffmann (2000) to explain the correlation. 
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Fig. 8. — Similar to Figure 7, but showing the effects of supermassive black holes. The dotted curves show results for star+halo models 
without black holes, for NFW (left panel) and Moore (right panel) models. The solid curves show results when central black holes are added. 
For the curves labeled "BH," the black hole masses are normalized by the empirical correlation between black hole mass and galaxy velocity 
dispersion (Ferrarese & Merritt 2000; Merritt & Ferrarese 2001). In the curves labeled "BHxA'^," the black holes are made systematically 
more massive by the factor N. (All galaxies have a = 1.0.) The heavy dashed curves again show the upper limits from six CLASS. 



5.3. Discussion 

When compared to the current hmits from CLASS, the 
star+halo models clearly predict too many detectable odd 
images. The discrepancy is not easily resolved by changing 
the density profile or invoking supermassive black holes. 
In other words, galaxies constructed from CDM mass dis- 
tributions have central densities that are too low. This 
conclusion is consistent with the result from Rusin & Ma 
(2001) that the lack of odd images requires steep den- 
sity profiles, p oc r~" with a > 1.8 at 90% confidence. 
Given that the measurements of odd images are only up- 
per limits, the constraints from odd images can only get 
stronger — and perhaps substantially, if more lenses turn 
out to have upper limits as strong as fodd < 0.0005 in 
B 0218+357. 

The odd image problem stands in stark contrast to most 
observational tests of CDM. Spiral galaxy dynamics arc 
said to imply that CDM halos are too concentrated to 
agree with observed galaxies (e.g., Moore 1994; de Blok et 
al. 2001; Salucci 2001; Weiner et al. 2001), although there 
is still substantial debate about whether beam smearing 
affects this conclusion for low surface brightness galax- 
ies (van den Bosch et al. 2000; van den Bosch & Swa- 
ters 2000). Independent of the dynamical arguments, the 
number and sizes of observed lenses lead to a similar con- 
clusion, as shown in §4. Recent interest in modifications 
to CDM, such as self-interacting dark matter (e.g., Spergel 
& Steinhardt 2000; see Wandelt et al. 2000 for a review), 
has therefore focused on making dark matter halos less 
concentrated. However, the lack of odd images implies 



that galaxies built from CDM halos arc not concentrated 
enough, to a significant degree. If self-interacting dark 
matter reduces densities, it would only exacerbate the odd 
image problem. 

It is clearly important to understand and resolve this 
paradox. Perhaps it is a case of comparing different sam- 
ples. Lensing intrinsically selects dense, massive galax- 
ies, so most lenses are elliptical galaxies; by contrast, the 
dynamical arguments against CDM halos come primarily 
from spiral galaxies, and in particular low surface bright- 
ness galaxies. However, the strongest limit on an odd im- 
age actually comes from a lens produced by a face-on spiral 
galaxy {fodd < 0.0005 for B 0218+357). 

More likely, it is a question of scales. Dynamical ob- 
servations and observed lensed images probe scales from 
several kiloparsecs down to ~0.5 kpc, while odd images 
probe scales more like tens of parsecs. The paradox could 
be resolved if halos have high densities on ~10 pc scales, 
low mean densities at ^1 kpc, and substantial dark matter 
halos beyond several kiloparsecs. The star+halo models 
do not show the required small-scale structure — but they 
are based on a questionable extrapolation of CDM profiles 
to scales much smaller than the resolution of numerical 
simulations. On such small scales, other effects may be 
important; while supermassive black holes appear not to 
suppress odd images at the required levels, self- interacting 
dark matter may be a mechanism to steepen the central 
cusp (Burkert 2000; Kochanek & White 2000; Moore et 
al. 2000), or simply to concentrate a lot of dark matter at 
very small radii (Ostriker 2000). Regardless of what the 
correct explanation turns out to be, it is clear that the odd 
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image problem provides a very interesting probe of galaxy 
mass distributions in the inner tens of parsecs. 

6. CONCLUSIONS 

Star+halo models of elliptical galaxies are a natural out- 
growth of the Cold Dark Matter paradigm that can repro- 
duce the quasi-isothermal mass distributions implied by 
stellar dynamics, X-ray halos, and gravitational lensing. 
When the stellar components arc fixed by observed galaxy 
populations, gravitational lens statistics can be used to 
constrain the dark matter components. The observed 
number and sizes of lenses place important upper limits 
on the amount of dark matter in the inner regions of ellip- 
tical galaxies: on average, dark matter can account for no 
more than about 33% of the total mass inside one effec- 
tive radius (Re), or about 40% of the mass inside 2i?e (95% 
confidence upper limits). Lensed images typically appear 
at a few effective radii, so the stellar and dark matter com- 
ponents must be comparably important for lensing.^ 

The dark matter limits are interesting when interpreted 
in the context of the CDM paradigm. Galaxies built from 
CDM mass distributions have significant amounts of cen- 
tral dark matter, so they predict that lenses should be 
more numerous and larger than observed. A high- density 
(Om = 1) CDM cosmology can therefore be ruled out at 
better than 95% confidence, while a low-density, flat cos- 
mology (r^M = 0.2, = 0.8) is at best marginally con- 
sistent with the lens data. By implying that CDM model 
galaxies have too much mass on kiloparsec scales, lensing 
independently supports the evidence against CDM from 
spiral galaxy dynamics (e.g., Moore 1994; de Blok et al. 
2001; Salucci 2001; Wcincr et al. 2001). The evidence from 
lensing is important because its only uncertainty is the de- 
composition into stellar and dark matter components, and 
lensing extends the argument to elliptical galaxies. These 
conclusions may be invalid if simple adiabatic contraction 
models do not apply to elliptical galaxies, but the mod- 
els appear to agree well with simulated galaxies even in 
merger scenarios (Gottbrath 2000). 

Lensing offers a unique additional insight into galaxy 
mass distributions. Star-|-halo models predict that a cen- 
tral or "odd" image should be detectable in more than 
30% of lenses, but such images are rarely observed. Equiv- 
alently, the upper limits on the fluxes of odd images for 
CLASS lenses lead to strong lower limits on the central 
densities of lens galaxies, and star-|-halo models fail to 
satisfy the limits. The failure is perhaps not surprising, 
because odd images are sensitive to the mass distribution 
on scales of tens of parsecs, where simple CDM models 
may break down. For example, supcrmassive black holes, 
which appear to be common in the centers of galaxies (e.g., 
Magorrian et al. 1998; Gebhardt et al. 2000a, 2000b; Fer- 
rarese & Merritt 2000; Merritt & Ferrarese 2001), might 
reconcile the models with the data by helping to suppress 
odd images (see Mao et al. 2001). However, star-|-halo 
models would need black holes that lie off the observed 
black hole-bulge mass correlation by more than a factor 
of 10 in order to suppress odd images at the required level. 
Alternatively, steeper central cusps could suppress odd im- 
ages, but only if mass cusps are considerable steeper than 

^ Dark matter mass fractions derived from individual lenses may differ 
averages and apply to the mass in spheres, while lens models give indi' 



luminosity cusps, only if steep cusps can survive merger 
events (which is unlikely if the progenitors contain black 
holes; Milosavljevic & Merritt 2001), and only at the ex- 
pense of aggravating the discrepancy in the number and 
sizes of lenses. 

The evidence from the number and sizes of lenses (and 
from spiral galaxy dynamics) implies, then, that CDM 
mass distributions have too much mass on kiloparsec 
scales; and the odd image problem indicates too little mass 
on tens of parsec scales. The problem may be the assump- 
tion that the dark matter particles are collisionless. Al- 
lowing dark matter self-interactions can lower the density 
on large scales (e.g., Spergel & Steinhardt 2000; Burk- 
ert 2000; Dave et al. 2001), although matching observa- 
tions may require substantial fine timing (e.g., Kochanek 
& White' 2000; Moore et al. 2000; Yoshida et al. 2000). As 
an intriguing corollary, Ostriker (2000) proposes that self- 
interactions could also increase the dark matter density 
on very small scales. Self-interacting dark matter models 
might therefore provide a way to reconcile cosmological 
mass models with real galaxies over a wide range of spa- 
tial scales, although the details remain to be worked out. 
In any case, it appears that lensing will be a very im- 
portant test of modifications to CDM through its ability 
to simultaneously test mass distributions on large scales 
(via observed images) and small scales (via limits on, and 
eventually detections of, odd images). These tests based 
on galaxy-scale lenses will complement constraints on self- 
interacting dark matter from giant arcs produced by clus- 
ter lenses (Meneghetti et al. 2001). 

It would be interesting to apply star-l-halo models to 
individual lenses, especially given the statistical limits im- 
plying that the stellar and dark matter components are 
of comparable importance in lensing. Comparing the 
halo/galaxy mass ratios inferred for individual lenses with 
the limits from statistics would be an important test of the 
models. Using two independent components would make 
it possible to to examine whether the stellar and halo dis- 
tributions have similar or different ellipticities and orien- 
tations, although the decomposition may not be unique. 
Finally, using two components with different ellipticities 
and/or orientations would test whether single-component 
lens models are sufficient or oversimplified. In particular, 
star-hhalo models might provide an internal reason why 
most lenses cannot be fit by a single ellipsoidal mass distri- 
bution (e.g., Keeton, Kochanek & Seljak 1997), although 
internal effects may often be smaller than external tidal 
perturbations from objects near the lens galaxy or along 
the line of sight (e.g., Hogg & Blandford 1994; Schechter 
et al. 1997). 

The constraints on dark matter from lens statistics can 
be strengthened with better data in at least three ways. 
First, the strongest constraints come from the distribution 
of image separations, where the model predictions are sen- 
sitive to the redshift distribution of the CLASS sample. 
The CLASS spectroscopic subsample (Marlow et al. 2000) 
provides a useful starting point, but the redshift distribu- 
tion of the full sample must be better constrained to make 
the lensing analysis truly robust. Second, many of the cur- 
rent constraints on odd images lie at the fodd — 0.01 level 

somewhat from the limits just quoted. The quoted limits are statistical 
ddual masses and involve the mass in cylinders. 
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(Rusin & Ma 2001). If the limits could be improved to the 
fodd — 0.003 level or better, they would make the mod- 
els much more sensitive to supermassive black holes, and 
lensing would become a powerful probe of black holes in 
distant galaxies out to rcdshift z ^ 1. Finally, this anal- 
ysis is based primarily on the sample of lenses from the 
CLASS survey, which is the largest existing lens survey 
but still has only 18 lenses. Larger surveys, in particular 
the Sloan Digital Sky Survey (SDSS; York et al. 2000), 
should increase the number of lenses by well over an order 
of magnitude. The SDSS lens sample will dramatically im- 
prove the constraints from lensing, provided that selection 
effects are well understood and that there is a subsample 
of radio loud lenses where useful limits on odd images can 



be obtained. 
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APPENDIX 

A. AN ANALYTIC SOLUTION OF ADIABATIC CONTRACTION 

Blumenthal ct al. (1986) give a simple analytic treatment of spherical adiabatic contraction that agrees remarkably 
well with more detailed numerical simulations. Let Mi(ri) be the initial mass profile as a function of the initial radius r^, 
while Mg{r) and Mh{r) are the final mass profiles of the galaxy and halo, respectively. In the Blumenthal et al. (1986) 
prescription, the three profiles are related by two equations, 

r [Mg{r) + Mh(r)] = n M^in) , (Al) 

Mhir) = {1 - fcooi)M,in) , (A2) 

where fcooi = Mg^tot/Mi^tot is the fraction of the system's mass contained in baryons that cool to form the galaxy. (There 
can be other baryons that remain hot and distributed throughout the halo, but they do not affect adiabatic contraction.) 
This adiabatic contraction formalism has often been applied to the problem of a disk galaxy in an NEW halo. Rix et 

al. (1997) have computed adiabatic contraction for elliptical galaxies numerically, but 1 find that with a Hernquist model 
(eq. 1 in §2.1) the problem can be solved analytically. Each initial radius maps to a unique final radius r given by the 
solution of the equation 

which is a cubic polynomial in r. Note that 1 have take the galaxy scale radius from eq. (1) and relabeled it as Sg. Also, 
= Mi{ri)/M2oo is the initial mass profile normalized by the virial mass (the mass inside the virial radius r2oo)- In 
the limit r ^ Sg, eq. (A3) has the simple asymptotic solution 

r= . (A4) 

fcool + (1 - fcool)mi{n) 

The full general solution can be also be found analytically, although it cannot be written quite so compactly. Following 
Abramowitz & Stegun (1981), solve a cubic equation of the form 

+ a^z^ + a\z -|- ao = (-^5) 

by defining 



fcool + (r + Sgf [(1 - fcooiV - n] TUiin) = , (A3) 



aia2-3ao (4 /a^x 
P= 5 ^. (A6) 

(AT) 

(p+v^^^y^\ (A8) 



S2 



{p-V^+?y^. (A9) 

zi = (si + 52) - y . (AlO) 

There are two other roots that may be real or complex, but because the ri r mapping under adiabatic contraction 
should be one-to-one, only the single real root is relevant. Once the cubic equation has been solved to map to r, eq. (Al) 
can be used to write the total mass profile as 

MUr)^Mgir)+M,ir) = ^-^MM)- (AH) 

This solution of adiabatic contraction by a Hernquist galaxy can be used for any form of the initial halo, by simply 
inserting the desired initial profile Mj(ri) into eqs. (A3) and (All). 



There is always a real solution of eq. (A5) at 
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B. THE MEAN MAGNIFICATION OF ODD IMAGES 

The mean magnification of odd images can be computed analytically, at least for simple lens models. Neglecting 
magnification bias, the mean magnification is defined to be 

{l^odd) = j^^^^ , (Bl) 

where the integrals extend over the multiply-imaged region of the source plane. Changing variables in the numerator to 
integrate in the image plane yields 

M = , (B2) 

where the integral in the numerator extends over the region in the image plane where odd images are found. The result 
is so simple because the figddi'^) factor in eq. (Bl) is exactly cancelled by the .Jaeobian of the coordinate transformation. 
Eq. (B2) says that the mean odd image magnification is simply the area where odd images occur in the image plane 
divided by the area of the multiply-imaged region of the source plane. This result is general and does not require specific 
symmetries in the lens model. It can be generalized to other types of images as well, provided that multiplicities are 
properly counted. 

Now focusing on spherical systems, {fiodd) — {fcr/uoutf' where Vcr is the radial critical curve and Uout is the boundary 
of the multiply- imaged region (see §3). These radii are found as follows. The critical radius is the solution of the equation 

= 1, (B3) 



dR 

where (pR is the lensing deflection. The boundary of the multiply-imaged region is then 

Uout = {cl>R-R)\ji^^^^. (B4) 

Consider two simple models. First, for a softened isothermal sphere with density poc {s'^ + r^)~^ with core radius s, 

As 

^^"''^ = (V4r, + .-3V.)^ ' (^'^ 



rs \rE J 2 yrEj \rEj 

where ve is the Einstein ring radius of the model when the core radius is zero. Second, for a simple power law density 
p oc r~" (with a > 1 to ensure that the projected mass distribution is a decreasing function of radius), 

(Mo..) = If 1)]' ;J^5«<2 (B7) 



Note that for a > 2 the model does not produce odd images (the density is singular, so the odd image theorem does not 
apply; see §5.1), so {podd) = 0. 
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